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Abstract 

The sensitivity of proton-induced J/^ and T production to the possible 
u/d asymmetry in the nucleon is studied. The ratio of the quarkonium 
production cross sections at large xp on hydrogen over deuterium targets, 
crip + p) / cr{p + d) , is shown to be sensitive to this asymmetry. Predictions of 
various theoretical models for this ratio are presented. 



The New Muon Collaboration (NMC) reported the result of their measurement of the 
Gottfried Sum ^ in muon deep inelastic scattering (DIS). The Gottfried Sum is given as 

/g(xi, X2; Q') = r dx{F2^P{x, Q') - F2^"(x, Q^))/x, (1) 

J Xl 

where F2^^ and F2^"' are the proton and neutron structure functions measured in muon DIS. 
Assuming charge symmetry for the parton distributions in the proton and neutron, F2 can 
be expressed in terms of the valence and sea quark distributions of the proton. The Gottfried 
Sum becomes 

lGixi,X2;Q'^) = - dx[uy{x,Q'^) - d^{x,Q'^)] + - I dx[u{x,Q'^) - d{x,Q'^)]. (2) 

O Jxi O J Xl 

The integral in the first term of (2) gives the difference between the number of the up and 
down valence quarks in the proton. The second term in (2) vanishes under the assumption 
of a M — J flavor-symmetric sea for the proton, and the Gottfried Sum Rule (GSR) 

lGiOA;Q') = l (3) 

is obtained. Based on their measurements of muon DIS on hydrogen and deuterium targets, 
the NMC has determined Jg(0.004, 0.8; = AGeV^) = 0.227 ± 0.007 Following an 
extrapolation to the unmeasured x-region, it was estimated /g(0, 1; = AGeV"^) = 0.240 ± 
0.016, significantly different from 1/3. 

Many theoretical models [3-13] have been proposed to account for the apparent violation 
of the GSR. These models in general fall into two categories. Models in the first category 
assume that the valence quark distributions in the proton are sufficiently singular at x < 
0.004 such that a large contribution to the Gottfried Sum occurs at this region not probed 
by the NMC experiment. Therefore, the GSR is not violated and the assumption oi u/d 
symmetry in the proton remains valid in these models. The other category of theoretical 
models [5-13] interpret the NMC result as evidence that the u and d distributions in the 
proton are different. Some empirical expressions for d{x) — u{x) have been proposed [^,^|13|. 
Several recent sets of parton distribution functions [|T^ , [l5[| explicitly allow u/d asymmetry 
to account for the NMC result. The origin of the enhancement of d over u in the proton 
has been attributed to pion cloud []5|,[7|J^,p!^ , diquark clustering in the nucleon [0, as well 
as Pauli-blo eking effect [|Tl[ . 



It has been proposed p, p!6| , p!7| that the Drell-Yan process provides an independent and 
sensitive test of the possible u/d asymmetry in the proton. Using a proton beam and 
restricting the kinema.tic reerime to forward xi? (xjt > 0.2V it is straiffhtforwaxd to show that, 



^ (Tdy{p + p) ^ ^ _ dp{x) - Up(x) 
(^Dvip + d) dp{x)+Up{x)' 

A similar expression with a reduced sensitivity to the u/d asymmetry can also be obtained 

for the ratio of Drell-Yan cross sections on neutron-rich targets versus isoscalar targets. In 



fact, the E772 Drell-Yan data [Il8| obtained with tungsten and isoscalar targets have been 



compared with predictions from various models. More recently, the NA51 experiment ||19 



reported 2cr£iy(p + p)/crDYip + d) = 0.91 ± 0.02 ± 0.02 measured at 450 GeV near xp ^ 
and X = 0.18. The NA51 result shows a significant asymmetry of u and d in the proton. 
Another Drell-Yan experiment covering a wider x range (0.05 < a; < 0.3) at 800 GeV has 
also been proposed pOl) . 

In addition to the DIS and the Drell-Yan processes, there are other interactions sensitive 
to the sea-quark distributions in the nucleon. Of particular interest is the hadronic produc- 
tion of J/ip and T, which can be detected in the same experiments designed to measure the 
Drell-Yan process. In this paper, the sensitivity of J/ip and T production to the sea quark 
distributions in the nucleon is studied. 

In contrast to the electromagnetic Drell-Yan process, hadronic J/ip and T production 
involves strong interactions. The two processes responsible for producing a pair of heavy 
quarks (cc and bb) are the qq annihilation and gg fusion. The cross sections for these QCD 
subpro cesses are pTI 



2, 37ra? 



a{qq ^ QQ; m') = i^{rn' + 2m^)A, 

2 2 \ 1 

a{gg - QQ; m^) = ^,{{m' + Am'ml + m%) H"^^) ' -^{^rn^ + 3K)A}, (5) 

where m is the invariant mass of the heavy quark pair {fQQ\ mq is the heavy quark mass 

and A = (m^ — 4m^mQ^)^/^. According to the QCD factorization theorem, the differential 

cross section for producing a QQ pair in hadronic interaction is 

da 2r 2\ 

d^ = (a:^2 + 4r^)V2^^^(^^'^^'^ )' (6) 

where Xi, X2 are the fractional momenta carried by the projectile (P) partons and target (T) 
partons, respectively, and = w? jS. HpT is the convolution of parton cross 

sections and the parton distribution functions in the projectile and target hadrons 

HpT{xi,X2]ni^) = Gp{xi)Gt{x2) (T{gg ^ QQ]m^) + 

{(lUxi)qkx2) + qp{xi)qTix2)}(r{qq^QQ]m^), (7) 



G{x),q{x), and q(x) signify the gluon, quark, and antiquark distribution functions respec- 
tively. 

To go from the production of QQ pair to the production of QQ bound states, we use the 
semi-local duality model ||22|. In this model, the QQ bound state cross section is obtained 
by integrating the free QQ cross section over r from the QQ threshold, ri = 2mQ/y^, to 
the open charm (or beauty) threshold, T2 = 2m£)(^B)/ \/s- Hence 

da/dxpiJ/ip, T) = F H 2TdTHpT{xi, X2; m'^)/{xF^ + 4x^)1/2, (8) 

Jti 

where F is the fraction of the QQ bound state cross section leading to J/ip{T) production. 

Despite its simplicity, the semi-local duality model is capable of describing many features 
of hadronic J/ip and T productions |23-p6l. In particular, the shape of da/dxp as well as 
the dependences of the cross sections on the beam energy and the incident-particle type, 
which are very sensitive to the relative contributions of gg fusion and qq annihilation, are 
well reproduced. The success of the semi-local duality model suggests that it gives a reliable 
description for the relative importance of these two processes. 

Figure shows the da/dxp data p6|j27[| for J/ip and T production using an 800 GeV 
proton beam. Results of the semi- local duality model calculations using Eq.(8) are shown as 
the solid curves. The quark masses nic and mb were set at 1.5 GeV and 4.5 GeV, respectively. 



in the calculation. The recent structure function set, DOl.l ||2^, was used. The shapes of 
the differential cross sections are well described by the calculations, and F is determined to 
be 0.17 and 0.034, respectively, for J/ip and T production. 

The contributions from gg fusion and qq annihilation to the J / if) and T production cross 
sections are shown as the dashed and dotted curves, respectively, in Figure |I|. For J /if) 
production at this beam energy, gg fusion is the dominant process. However, at xp > 0.6, 
qq annihilation process starts to dominate. This reflects the fact that the gluon distribution 
drops more rapidly at large x than does the quark distribution. For T production. Figure |l] 
shows that qq annihilation contributes ~35% of the cross sections near xp = 0, and it is the 
dominant contribution at large xp. Since T is roughly three times more massive than J/ip, 
production of T is more sensitive to structure functions at large x, where gluon ceases to 
dominate the parton densities. This accounts for the relative importance of qq annihilation 
for T production. 

For J/ and T production in p — p interactions, the cross section from uu and dd anni- 
hila.tinn. rr(n-\-n). is nrnnnrtinnal tn 



Up{Xi)Up{x2) + Up{Xi)Up{x2) + dp{xi)dp{x2) + dp{xi)dp{x2). (9) 
For p — d interactions, the corresponding cross section, a{p + d), is 

Up{Xi){Up{x2) + Unix2)) + Up{Xi) {Up{x2) + M„(a;2)) 

+dp{xi){dp{x2) + dn{x2)) + dp{xi){dp{x2) + dn{x2)). 

In Eq.(lO) the deuteron parton distributions are taken as the sum of the proton (p) and 
neutron (n) distributions. Charge symmetry requires Up{x) = dn{x),Up{x) 
Expression (10) now becomes 

Up{Xi){Up{x2) + dp{x2)) + Up{Xi){Up{x2) + dp{x2)) 
+dp{xi){dp{x2) + Up{x2)) + dp{xi){dp{x2) + Up{x2)) 

For a flavor-symmetric sea, Up = dp, Eqs. (9) and (11) give a{p + d) = 2a{p + p). This 
relation is also valid for contributions from the ss annihilation and gg fusion. Therefore, the 
ratio R{xf), defined as 

nr^ ^ _^ da/dxF {p + p^ J/i^jT)) , . 

~ da/dxF{p + d^ J/^(T)) 

would be equal to 1 for all models which assume Up = dp. On the other hand, R{xf) could 
deviate significantly from 1 if -Up 7^ dp, especially at large xp where qq annihilation has a 
dominant contribution. 

Figure ^ shows the predictions of R{xf) for J/tp and T production with 800 GeV proton 
beam. Five different structure function sets have been used in the calculations. The solid 
curves correspond to the u/d symmetric DOl.l structure functions. As discussed earlier, 
R{xf) is identically one in this case. The dashed and dotted curves are results for the 
u/d asymmetric structure function sets MRSD-' and CTEQ2pM, respectively. These two 
structure function sets were obtained from recent global fits to Drell-Yan and DIS data 
including the NMC result |n,T^J2^. Finally, the calculations using the parameterization of 



d{x) — u{x) given by Ellis and Stirling and Eichten et al. ||T3| are also shown in Figure |[ 
The parameterizations for d{x) — u{x) were given at a fixed value of at 4 GeV^. The Q^- 
dependence was taken into account in the calculation by using the Altarelli-Parisi evolution 
equation for flavor non-singlet structure functions [Q. 

Figure ^ shows that the measurement of R{xf) for T production at 800 GeV provides 
a sensitive test of models which have different descriptions on the u/d symmetry in the 
proton. In comparison, J/ip production is much less sensitive to the sea quark distributions 



due to the dominance of gg fusion at this energy. However, at lower beam energies, J /if) 
production is more sensitive to the sea quark distributions as it probes a larger x region 
where qq annihilation starts to be important. This is illustrated in Figure |^ which shows the 
predictions of R{xf) for J/tp production using 450 GeV proton beam. The recent result from 
NA51 is also compared with the predictions. While the predictions are consistent with 
the data, it is important to extend the measurements to larger xp values. A measurement 
of R{xf) for J/tp production at even lower beam energies would also be of interest. 

It should be mentioned that a more elaborate scheme for calculating hadronic J/\E' and 
T production, the so-called color-singlet model, has also been proposed [|3ll , p2| . This model 
includes relevant gg, gq, gq, and qq processes up to order of a^, and is capable of predicting 
the absolute production cross sections. It would be very interesting to compare the pre- 
dictions of the semi-local duality model and the color-singlet model with the experimental 
data. 

In conclusion, proton-induced J ftp and T production offers an independent means to 
examine the possible u/d asymmetry in the proton. Measurement of the cross section ratios 
on hydrogen and deuterium targets in the forward xp region would provide a sensitive test 
of various theoretical models on the sea quark distributions of the proton. 



REFERENCES 



[1] P. Amaudruz et ai, Phys. Rev. Lett. 66 (1991) 2712. 
[2] K. Gottfried, Phys. Rev. Lett. 18 (1967) 1174. 

[3] A. D. Martin, R. G. Roberts, and W. J. Stirling, Phys. Lett. 252B (1990) 653. 
[4] P. N. Harriman, A. D. Martin, W. J. Stirhng, and R. G. Roberts, Phys. Rev. D42 (1990) 
798. 

[5] E. M. Henley and G. A. Miller, Phys. Lett. 251B (1990) 453. 
[6] S. D. Elhs and W. J. Stirhng, Phys. Lett. 256B (1991) 258. 

[7] S. Kumano, Phys. Rev. D43 (1991) 59, S. Kumano and J. T. Londergan, Phys. Rev. 
D44 (1991) 717. 

[8] A. Signal, A. W. Schreiber, and A. W. Thomas, Mod. Phys. Lett. A6 (1991) 271. 

[9] G. Preparata, P. G. Ratcliffe, and J. SofFe, Phys. Rev. Lett. 66 (1991) 687. 
[10] M. Anselmino and E. Predazzi, Phys. Lett. 254B (1991) 203. 
[11] R. D. Field and R. P. Feynman, Phys. Rev. D15 (1977) 2590. 
[12] H-Y. P. Hwang and J. Speth, Phys. Rev. D46 (1992) 1198. 
[13] E. J. Eichten, I. Hinchhffe, and C. Quigg, Phys. Rev. D45 (1992) 2269. 
[14] A. D. Martin, R. G. Roberts, and W. J. Stirling, Phys. Lett. 306B (1993) 145. 
[15] J. Botts et al, Phys. Lett. 304B (1993) 159. 

[16] S. Kumano and J. T. Londergan et ai, Phys. Rev. D46 (1992) 457. 
[17] W-Y. P. Hwang, G. T. Garvey, J. M. Moss, and J. C. Peng, Phys. Rev. D47 (1993) 
2697. 

[18] P. L. McGaughey et al, Phys. Rev. Lett. 69 (1991) 1726. 
[19] A. Baldit et al, to be published (1994). 

[20] G. T. Garvey et al, Fermilab Proposal E866 (unpublished) (1992). 
[21] B. L. Combridge, Nucl. Phys. B151 (1979) 429. 
[22] H. Fritzsch, Phys. Lett. 67B (1977) 217. 

[23] V. Barger, W. Y. Keung, and R. J. N. PhiUips, Z. Phys. C6 (1980) 169. 

[24] S. Falciano et al, Phys. Lett. 158B (1985) 92. 

[25] R. Vogt, S. J. Brodsky, and P. Hoyer, Nucl. Phys. B360 (1991) 67. 

[26] M. S. Kowitt et al, Phys. Rev. Lett. 72 (1994) 1318. 

[27] G. Moreno et al, Phys. Rev. D43 (1991) 2815. 

[28] J. F. Owens, Phys. Lett. 266B (1991) 126. 



[29] H. Plothow-Besch, "PDFLIB : Nucleon, Pion and Photon Parton Density Functions and 

as Calculations," CERN-PRE, 1994.01.11 (1994). 
[30] G. Altarelli and G. Parisi, Nucl. Phys. B126 (1991) 298. 
[31] R. Baier and R. Ruckl, Z. Phys. C19 (1983) 251. 

[32] R. Gastmans, W. Troost, and T. T. Wu, Nucl. Phys. B291 (1987) 731. 



FIGURES 



10' 



J2 
X 



10 



b 



1 



p + Cu ^ J/Y' 

800 GeV 








III 


1 1 1 1 1 1 1 1 1 1 

p + Cu ^ T 


II- 

b) : 




800 GeV 






^^^^^^^^^^ 










III 


1 1 1 1 1 1 1 1 1 1 1 1 1 


N _ 
N — 

1 



-0.6 



-0.4 



-0.2 







0.2 



0.4 



0.6 



FIG. 1. J/V-' and T production cross sections inp+Cu interactions from References [26], [27]. 
The solid curves are calculations using the semi- local duality model (Eq.(8)) and DO 1.1 structure 
functions [28]. The contributions from the gg and qq annihilation processes are shown as dashed 
and dotted curves respectively. 
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FIG. 2. Predictions of R{xf) for J/V' and T production at 800 GeV using various structure 
function sets in the semi-local duality model. 
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FIG. 3. Predictions of R{xp) for J/V' production at 450 GcV using various structure function 
sets in the semi-local duality model. The data from NA51 [19] is also shown. 



